Doping effect on thermoelectric properties of M0S2 
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We systematically study the potential of layered M0S2 as a promising thermoelectric material for 
the first time, based on Boltzmann transport theory and first-principles calculations. Anisotropic 
thermoelectric transport property are investigated and compared with experimental data. A good 
agreement between our calculated temperature-dependent thermopower and available experimen- 
tal data is obtained. We then discuss the temperature and doping level dependent thermopower, 
electrical conductivity, power factor (PF) and ultimately figure of merit (ZT) coefficient along in- 
plane and cross-plane directions. M0S2 has a vanishingly small anisotropy of thermopower but a 
big anisotropy of electrical and thermal transport properties in the optimal doping region. A clear 
preference for thermoelectric transport under electron doping, as well as for the cross-plane trans- 
port direction is found. Our study demonstrates that it is essential to improve the thermoelectric 
performance of M0S2 by reducing thermal conductivity through restacking exfoliated layers, or by 
electron doping through substituting at S site. 



I. INTRODUCTION 

Thermoelectrics play a key role for power genera- 
tion and refrigeration^^. Chalcogenide composite ma- 
terials, such as Bi2Te3, PbTe and others^^^, have at- 
tracted much research interest for decades because of 
their high performance on thermoelectric conversion. 
Transition- metal dichalcogenide M0S2, one prototype 
material in chalcogenide family, has distinctive electronic, 
mechanical, catalytic and tribological properties^ 
Recent research on its optical properties^! and lattice 
dynamics^^^^^ has aroused renewed interest but its po- 
tential thermoelectric capability has barely been stud- 
ied. Nevertheless, the lowest thermal conductivity ob- 
tained experimentally in a MoS2-related structure^^, to- 
gether with an unusually large thermopower^^^^ found 
in M0S2, may render it a potential candidate for thermo- 
electric application. 

Unfortunately, available experimental data concerning 
the thermoelectric transport properties of M0S2 are quite 
scarce and fragmentary in literature, and no theoreti- 
cal studies have thus far been found. The most signifi- 
cant experimental work was performed by Mansfield and 
Salam^^ and Thakurta et. al.^^, who studied tempera- 
ture dependence of the electrical transport properties in- 
cluding thermopower S and electrical conductivity a at 
only three samples with low dopings, but failed to inves- 
tigate on the thermal-related properties. Kim et. al.^^ 
merely worked on the thermal conductivity Based on 
these incomplete experimental data, it is difficult to cal- 
culate the figure of merit ZT which is crucial to evaluate 
thermoelectric conversion capability. Moreover, available 
extrinsic carrier concentrations in M0S2 single crystals 
were quite low. It remains a big challenge to obtain a 
wide doping region in experiment but nonetheless is es- 
sential to tuning the carrier concentration appropriate 
for an optimum ZT. 

In this letter, we study theoretically the thermoelec- 



tric transport properties of M0S2 over a wide range of 
doping level (from 10^^ to 10^^ cm~^) and focus on 
their directional anisotropics. We have the following 
findings: (1) Thermopower is attainable more than 200 
/iV/K over a wide range of dopings and directional ther- 
mopower anisotropy vanishes between in-plane and cross- 
plane thermopowers for doping level above 10^^ cm~^. 

(2) Anisotropic electronic scattering time exists between 
in-plane and cross-plane directions, which accounts for 
two orders of magnitude difference between (Jxx and Gzz • 

(3) A clear preference for thermoelectric transport under 
electron doping, as well as for the cross-plane transport 
direction is found. (4) Heavy doping is essential for an 
optimal thermoelectric performance, the figure of merit 
ZT reaches 0.006 at TOOK and 4x10^^ cm-^. (5) ZT 
can be further optimized by reducing cross-plane ther- 
mal conductivity through random re-stacking of exfoli- 
ated layers. 



II. METHODOLOGY 

The band structure of M0S2 is calculated by us- 
ing the general potential linearized augmented plane- 
wave (LAP W) method as implemented in the WIEN2K 
package^^. The electronic exchange-correlation is de- 
scribed within the generalized gradient approximation 
(GGA) of Perdew-Burke-Ernzerhof(PBE) fiavor^^. We 
use 5000 k points in the full Brillouin zone(BZ) to achieve 
a self-consistency accuracy better than 1 meV/atom. 
M0S2 has Pea/mmc space group symmetry and consists 
of a hexagonal plane of Mo atoms sandwiched by two 
hexagonal planes of S atoms. The unit cell contains 
two alternating and weakly van-der-Waals-bonded lay- 
ers with an AB stacking along c axis. The experimental 
lattice parameters'^ (a = 3.16A, c = 12. 295 A) are used 
here. 

An improved band gap is essential for calculating 
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FIG. 1. (Color online) EV-GGA electronic band structure 
of M0S2 along high- symmetry lines (a) F-M-K-F in plane 
and (b) F-A across plane in the hexagonal Brillouin zone 
(BZ). The valence-band edge is set as zero and marked with a 
dashed line. Electronic density of states in (c) shows a com- 
paratively higher value at the conduction band edge than at 
the valence band edge. 



transport quantities, so the Engel-Vosko GGA (EV- 
GGA) formalism^^ is applied to calculate the band gap 
more accurately. In Fig. [H we present our calculated 
band structure and density of states of M0S2. The in- 
plane and cross-plane band structure are considered sep- 
arately in Fig. [TJa, b). In Fig. HJa), an indirect gap of 
1.04 eV is obtained between top of valence band at F and 
bottom of conduction band at one k point from K to F. 
A similar band structure calculated from standard PBE- 
GGA formalism gives an indirect gap of 0.84 eV, which 
agrees with reported valued. Compared with ~ 1.20 
eV from experiments^, it is clear that EV-GGA does 
improve the band gap calculation upon PBE-GGA. In 
contrast to the pronounced dispersive in-plane electronic 
bands shown in Fig.lTJa), the cross-plane bands shown in 
Fig. [TJb) are quite flat, showing a very weak cross-plane 
bonding due to Van der Waals interaction but a possi- 
bly better TE performance than in-plane bands^^. The 
structural anisotropy induces an anisotropy between in- 
plane and cross-plane band gaps. The size of calculated 
cross-plane gap is found to be 2.2 eV. 



III. RESULTS AND DISCUSSION 

Besides the band gap anisotropy, a strong asymmet- 
ric feature of band structure between valence and con- 
duction bands implies that the thermoelectric proper- 
ties of p-type M0S2 would be very different from that 
of n-type. The heavy and doubly degenerate bands near 
the conduction-band minimum suggest that the n-type 
M0S2 would have better thermoelectric performance. In 
Fig. [U^c), total density of states (DOS) also shows this 
preference, there is a comparatively higher DOS very 
close to the conduction band edge than that near the 
valence band edge. Considering more experimental data 
are found for p-type in literature^^^^, we focus on hole- 
doped M0S2 in this study. 
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FIG. 2. (Color online) (a) Temperature dependence of cal- 
culated in-plane thermopower Sxx of M0S2, compared with 
experimental data by Mansfield and Salam^^ at two doping 
levels p = 1.6 X 10^^ (filled triangle) and 3.4 x 10^^ (empty 
circle) holes per cm^. Doping level dependence of (b) in- 
plane thermopower Sxcc(p, T), (c) ratio of cross-plane Szz{p, 
T) over in-plane Sxcc(p, T) and (d) in-plane Sxx{n,T) at dif- 
ferent temperatures. The temperature ranges from lOOK to 
TOOK for some practical reason. Hole and electron doping are 
respectively used in (a-c) and (d). Experimental data with p 
= 7.6 X 10^^ and 1.0 x 10^'^ cm"^ at 200K by Mansfield and 
Salam^^ is respectively marked by filled triangle and square 
in (b). 



We calculate transport properties based on Boltzmann 
transport theory applied to the band structure. In 
the following part, we discuss the dependence of trans- 
port functions including S, a, PF, hc and ultimately ZT 
on temperature and doping level along two directions 
perpendicular to each other. The electron scattering 
time is assumed to be independent of energy due to its 
good description of S(T) in a number of thermoelectric 
material^^ i^^i^^ . In this sense, no adjustable parameters 
are needed to calculate those thermoelectric parameters. 
The integration is done within the BOLTZTRAP trans- 
port code^. A very dense meshes up to 18000 k points 
in the BZ is used. 

In-plane and cross-plane thermopower are firstly stud- 
ied. We initially compare our calculated S with experi- 
mental data from Mansfield and Salam^^ in Fig. [2fa,b). 
Samples in their experiment were p-type from nature 
with hole carrier concentration as low as 10^^ ~ 10^^ 
cm~^. A very good agreement is obtained. Both our 
calculation and available experimental data show a value 
of in-plane thermopower S higher than 400 /iV/K and 
S decreases with increasing temperature as shown in 
Fig.[2ja). Then we extend our discussion to high doping 
of 10^'^ ~ 10^^ cm"^ (corresponding to p = 10"^ ~ 10"^ 
holes per unit cell in our case) to optimize thermoelectric 
properties, as predicted from theory and as observed 
in many materials^^^^ including M0S2 that we discuss 
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in this letter. Therefore from now on, we put more 
attention on this high doping region. From Fig. Efb), 
the thermopower, though decreasing with doping level, 
takes a value of at least 200 /iV/K (a ZT~2.4 expected 
if Wiedemann- Franz law is assumed.) and increases with 
increasing temperature. To see a possible anisotropy be- 
tween in-plane and cross-plane S, we show the ratio of 
Szz over Sxx in Fig. Efc). Contrary to a relatively high 
anisotropy of thermopower at low doping level, little dif- 
ference is found at high level of doping. Not only the 
magnitude, but also the dependence on hole doping level 
and temperature are similar between S^^^^ and S^^. This 
is different from layered conductive thermoelectric ox- 
ides such as Naa^Co02 and Ca3Co409_(5^^'^^. To confirm 
that n-type M0S2 may be better, we briefly compare the 
data of n-type in Fig. [2jd) with p-type in Fig. EJb). Ex- 
pectedly, bigger value of thermopower of n-type M0S2 is 
found. 

Based on the energy-independent scattering time ap- 
proximation, it is quite straightforward to get doping de- 
pendence of (j/r at various temperature from the elec- 
tronic band structure, upon which we calculate a if r is 
known. In Fig.[3l^a), we show calculated Oxxj^xx depend- 
ing on hole doping and temperature. We find an almost 
temperature independence but an approximately linear 
doping level dependence of (Jxx/'Txxi namely (Jxx/'Txx ^ 
T^ p. For a quadratic band dispersion in an electron- 
phonon approximation, a ~ ps T~^ stands, this re- 
sults in Txx ^ T~^p~3, in consistent with the analyti- 
cal treatment of carriers scattered by lattice vibrations 
in a semiconductor— 1^. To calculate (j^xi we need some 
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FIG. 3. (Color online) Doping level dependence of (a) in-plane 
Cccx/tccx, (b) in-plane electrical conductivity ctxx, (c) ratio of 
cross-plane (Tzz/tzz over in-plane (Txx/txx and (d) cross-plane 
<Jzz^ at different temperatures. Isotropic and anisotropic elec- 
tronic scattering time r are respectively assumed in (c), for 
instance, Tzz/txx = 1 and 0.024 with the anisotropic one fit- 
ting well the experimental data in blue square from Thakurta 
et. aL— and thereby used to derive Gzz in (d). 



experimental input besides rxx{T,p). Here we use axx 
= 0.16 Q-^cm-^ at 100 K and 1.4x10^^ cm-^ dop- 
ing level from the experimental data by Thakurta et. 
al.^^. Our calculated Cxx/^xx combined with the exper- 
imental data yields Txx = 2.72 x 10~^^ s at lOOK and 
1.4x10^^ cm~^ doping level. Plugging it into our cal- 
culated cr/r, we show (Jxx{T^p) in Fig. [3^b), and power 
factor crxxSl^{T,p) in Fig.^a). 

Electrical conductivity along c axis is also calculated 
and a strong anisotropy is found between the in-plane and 
cross-plane electronic scattering time r. When isotropic 
scatter time r is assumed, namely, Tzz/txx = 1, we obtain 
^zz/^xx close to unity in Fig. [3](c), which is against the 
reported result—. To fit (Jxx/(^zz of two orders of magni- 
tude in experiment, we use anisotropic Tzz/txx = 0.024, 
and therefore obtain Tzz = 6.52 x 10~^^ s at lOOK and 
1.4x10^^ cm~^. This suggests that strong anisotropic 
scattering effects should play a role in this system. Fi- 
nally, we are able to calculate azz{T^p) and show it in 
Fig.Hd). 

With thermopower and electrical conductivity, we are 
able to evaluate power factor. For an optimized TE per- 
formance, peak value of PF and the corresponding dop- 
ing level and temperature are more concerned here. The 
power factors along two perpendicular directions both 
have peak values spanning in a wide doping range from 
10^^ to 10^^ cm~^. From Fig.HJ the value of peak power 
factor is nearly constant, while its temperature increases 
with the hole doping level. Nearly 50-fold difference is 
found between in-plane and cross-plane power factors, 
e.g., the in-plane and cross-plane peaks aS'^ at 700 K are 
respectively 4.1 x lO""^ W/m/K^ and l.Ox 10"^ W/m/K^. 
This has something to do with the anisotropic electronic 
scattering time r. Unfortunately, both values are rela- 
tively lower than that expected for a good thermoelectric 
material. An improvement on the inter-layer electrical 
conductivity seems very essential to improve PF. 

To optimize ZT value, the thermal conductivity a^, 
including electronic and lattice hci^ is also essential. 
Based on the scattering time approximation previously 
discussed, we calculate electronic thermal conductivity 
(shown in Fig. [5l(a,b)). At least two orders of magnitude 
difference between and hc^^ are found, which is in 
consistent with electrical conductivity case and fits well 
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FIG. 4. (Color online) Doping level dependence of (a) power 
factor axx Sxx and (b) (JzzSzz of p-type M0S2 at different tem- 
peratures. 
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FIG. 5. (Color online) (a) In-plane electronic thermal con- 
ductivity and (b) cross-plane as function of doping 
and temperature, (c-d) Temperature dependence of thermal 
conductivity Kzz from Kim et. al.^^ in experiment and Kxx of 
p-type M0S2. Kxx ^ SOOkzz is assumed in (d). The experi- 
mental data in (c) are presented by open and solid circles for 
pristine and exfoliated & restacked compound, respectively, 
with the pristine data fitted by Kzz = 183.103/T + 0.412671 
in dashed line. 
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FIG. 6. (Color online) Doping level dependent (a) in-plane 
and (b) cross-plane figure of merit ZT of p-type M0S2 at dif- 
ferent temperatures. 

directions, which is shown in Fig. [6fa,b). Optimum ZT 
value increases with increasing temperature, so does the 
corresponding optimum doping level. In-plane figure of 
merit falls behind cross-plane one by at least one order 
of magnitude, implying a better performance achieved 
along c-axis. However, ZT value is far below unity. To en- 
hance the ZT along c axis, a heavier electronic/hole dop- 
ing is desired and therefore upshifts its optimum working 
temperature. Besides doping, random stacking of two- 
dimensional crystalline layer to reduce the cross-plane 
thermal conductivity is also favorable. Kim et. al.^^ ex- 
foliated M0S2 films and then randomly restacked them, 
achieving a very low lattice k, (less than 0.7 W K~^ m~^ 
compared with the pristine crystal structure). But still 
there is some space to further reduce according to 
Chiritescu et. al.^^ where 30-fold reduction of thermal 
conductivity from crystal structure was obtained. 



IV. CONCLUSION 



with the Wiedemann-Franz law. For the lattice thermal 
conductivity, it is not possible to calculate f<ii from elec- 
tronic band structure and here we have to get it from 
experiment. We only find in literature the cross-plane 
Kzz from Kim et. al.^^ and show it in Fig. [5jc). The lat- 
tice thermal conductivity dominates, for instance, total 
hZzz = 1.00 W K~^ m~^ at 300K, which is two orders of 
magnitude bigger than let alone at lower tempera- 
ture. Obviously, hci is approximately equal to k,. We fit 
the experimental data for pristine M0S2 in Fig. [5jc) by 
using hi = 183.103/T + 0.412671. Noticing that Umklapp 
process has a 1/T dependence and thermal conductivity 
usually gets softening with increasing temperature. No 
experimental data exist for in-plane lattice a^, inspired 
by the graphite data^^, we assumed a 300 times bigger 
of Kxx than Hizz- 



To improve the thermoelectric conversion of M0S2, we 
suggest to make a heavy doping. However, it seem still 
insufficient since power factor is behind some standard 
value (10~^ W/m/K^). The problem lies in the poor 
electrical conductivity possible due to a slightly big elec- 
tronegativity difference between Mo and S elements. We 
suggest to increase a by substituting Se or Te for S. This 
has multiple advantages, e.g. to reduce the electron scat- 
tering induced by the electronegativity difference, to slow 
down the phonon propagation by using heavier atom, and 
in the meantime to realize the doping on M0S2. More en- 
couragingly, this may even enhance the thermopower— . 
Simply reducing further the thermal conductivity, espe- 
cially the lattice thermal conductivity without affecting 
power factor of M0S2, seems also possible as demon- 
strated on WSe2^^ and future experiment may be needed 
to test this possibility. We also have to mention that 
our discussion is limited below BOOK to avoid bipolar 
excitation^^, which is expected to increase thermal con- 
ductivity despite enhancing the electrical conductivity. 
How this bipolar excitation affects thermoelectrics above 
800K also requires a further study. 



All the data obtained above allow us to calculate ZT 
as function of temperature and hole doping along two 
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